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Description 

[0001] This invention relates to an apparatus for extracting features of skin pattern images and an image processor 
for smoothing and enhancing features of the skin pattern images making use of the extracted features. More particularly, 
5 the invention relates to an apparatus for extracting such features as ridge line directions, ridge line pitches or their 
combination of skin pattern images for fingerprint verification, fingerprint classification or palm print verification, and 
an image processor for processing skin pattern images for enhancing ridge lines by eliminating noises of the images 
making use of the extracted features. 

[0002] As a prior art, there is a Japanese patent application entitled "An Apparatus for Striped Pattern Extraction" 
10 and laid open as a Provisional Publication No. 77138/76. In the prior art, which is hereafter called a first prior art, an 
apparatus is disclosed for eliminating noises from such striped images as fingerprints or tree rings and extracting both 
stripe patterns and directional patterns of the striped images. 

[0003] FIG. 14 shows principle of the apparatus. An image in a subregion 100 with its center at a small square 101 
is represented by a set of pixel data. Each small square in the subregion 100 represents a pixel, and a numeral 1 or 
15 a blank in a small square means corresponding pixel value 1 (black) orO (white). Filter data sets 102, 103, 104, 105, 
and 1 06 for masking are prepared corresponding to no direction dO and four directions tfl , d2, d3, and ctt respectively. 
By calculating their sum of products with the pixel data of the image exampled in FIG. 14, values +5, 0, +11, +3 and 
-7 are obtained corresponding to each direction dO to d4 respectively. 

[0004] From these values, ridge direction at pixel 1 01 of the image is determined to the direction <J2 which gives the 
20 maximum value +11 and stripe pattern value at the pixel 101 is determined to the maximum value +11 . 

[0005] Thus, stripe patterns and directional patterns are extracted as features of striped images in the first prior art. 
[0006] Smoothing of extracted features is sometimes necessary for decreasing noise influence on the features. For 
the purpose, a smoothing method is disclosed in a Japanese patent application entitled "Method of Smoothing a Ridge 
Direction Pattern and an Equipment used therefor" and laid open as a Provisional Publication No. 18195/93. This 
25 smoothing method, hereafter called a second prior art, is based on a minimum energy principle. 

[0007] From ridge direction data extracted with their reliability coefficients, a smoothed ridge direction pattern is 
determined in the prior art, wherein it is assumed that direction data extracted with the higher reliability need the more 
energy to be smoothed into a direction pattern and that it needs the more energy to smooth direction data into a direction 
pattern the more different from neighboring direction patterns. Therefore, a direction pattern which gives a minimum 
30 value of an evaluation function, linear summation of those energies, is determined as the smoothed direction pattern 
of a subregion. 

[0008] However, pixel values of a striped image vary periodicaly in a direction perpendicular to the stripes. Therefore, 
when stripe phase of a striped image is nearly orthogonal to the stripe phase of the filter data for masking used in the 
first prior art, the sum of products results nearly 0 even when both directions coincide with each other. This means that 
35 the extracted direction patterns and, consequently, the extracted stripe patterns become unreliable influenced by noise 
data, making precise ridge line pitch extracion difficult. 

[0009] In the second prior art, smoothing of ridge direction patterns is performed based on extracted ridge directions. 
So, when extracted ridge direction data are not reliable, precise smoothing of ridge direction patterns and, consequently, 
precise smoothing of ridge line pitches can not be performed. 

40 [0010] Another prior art document can be found in IEE PROC.-VIS. IMAGE SIGNAL PROCESS., vol. 141 , no. 2. 
April 1994, pages 87*94, XP000610431 B. G. SHERLOCK, D M. MONRO ETK. MILLARD: "Fingerprint enhancement 
by directional Fourier filtering". In this document, a new method of enhancing fingerprint images is described, based 
upon rionstationary directional Fourier domain filleting. Fingerprints are first smoothed using a directional filter whose 
orientation is everywhere matched to the local ridge orientation. Thresholding then yields the enhanced image. 

45 [001 1] It is the object of the present invention to overcome the drawbacks of the prior art documents. 
[0012] The object is solved by the apparatus according to independent claims 1 and 2. 

[0013] Further advantageous features of the present invention will become apparent by the dependent claims, the 
description and the accompanying drawings. 

[0014] The invention provides an apparatus for extracting features such as ridge directions, ridge line pitches or their 
50 combination from skin pattern images with high reliability regardless of the stripe phase difference, and a skin pattern 
image processor for precisely enhancing ridge line patterns by eliminating noises from the skin pattern images. 
[0015] Furthermore, a plurality of filters consisting of two-dimensional weight coefficient matrix are prepared in the 
present invention according to kinds of features to be extracted. 

[0016] A skin pattern image is filtered with each by each of the plurality of filters by convolution calculation on a real 
55 plane or by product calculation of each frequency component on a Fourier transformed plane. 

[0017] Contrast or image intensity, that is, a square sum of pixel values of a subregion of each filtered image thus 
calculated is compared with each other and a feature parameter represented by a filter which gives a largest image 
intensity is determined as the feature to be extracted of the subregion of the skin pattern image. 
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[0018] And from pixel values of the image data of subregions filtered by a filter giving the feature of the subregion, 
image data of the subregion with enhanced feature can be obtained at the same time. 

[0019] The foregoing, further aspects, features, and advantages of this invention will become apparent from a con- 
sideration of the following description, the appended claims, and the accompanying drawings in which the same nu- 
5 merals indicate the same or the corresponding parts. 

[0020] FIG. 1 is a block diagram illustrating a first embodiment of the invention. 

[0021] FIG. 2 illustrates an example of a filter designed on a Fourier transformed plane. 

[0022] FIG. 3 shows frequency bands of five filters prepared for extracting ridge line pitches. 

[0023] FIG. 4 shows shading images of the five filters of FIG. 3. 
10 [0024] FIG. 5 illustrates an example of arrangement of subregions in a skin pattern image. 

[0025] FIG. 6 is a block diagram of a filtering means 24 of a second embodiment of the invention. 

[0026] FIG. 7 is a block diagram of a feature extracting means 30 of a third embodiment of the invention. 

[0027] FIG. 8 illustrates a subregion C(m, n) and neighboring subregions C(m' t n') around the subregion C(m, n). 

[0028] FIG. 9 is a flowchart illustrating operation of the feature extracting means 30 of FIG. 7. 
15 [0029] FIG. 10 is a block diagram illustrating a fourth embodiment of the invention. 

[0030] FIG. 11 is a block diagram illustrating a fifth embodiment of the invention. 

[0031] FIG. 12 shows shading images of examples of filters prepared in a first filter file 811 of FIG. 11 . 

[0032] FIG. 13 shows shading images of examples of filters prepared in a second filter file 812 of FIG. 11 . 

[0033] FIG. 14 is a diagram illustrating principle of a prior art. 
20 [0034] Now, embodiments of the present invention will be described in connection with the drawings. 

[0035] Referring to FIG. 1 , a first embodiment of the invention comprises a image memory 11 for storing image data 

of such skin patterns as fingerprint pattern images, a filter file 1 5 for storing filters prepared for filtering the image data 

read out from the image memory 1 1 , filtering means 12 for filtering the image data by convolution calculation each by 

each with filters stored in the filter file 15, image intensity calculation means 13 for calculating image intensities of 
25 subregions of each of the filtered data, and feature extracting means 14 for determining features of predetermined 

subregions of the image data, in accordance with the results calculated in the image intensity calculation means 13. 

[0036] The image memory 11 is composed of a hard disk memory or a DRAM for storing digital image data of skin 

pattern images as fingerprint pattern images or palm print pattern images collected by an input device as a CCD camera 

or an image scanner. 

30 [0037] In the filter file 15, a plurality of filters used by the filtering means 12 are prepared according to a kind of 
features to be extracted in a DRAM for example. 

[0038] In the following descriptions, a set of image data stored in the image memory 11 is assumed to be data of 
512 x 512 pixels with resolution of 20 pixels per mm, and is expressed by a function gfx.y) representing pixel values 
at coordinates (x,y), where x = 0, 1,2,511 andy=0, 1,2,..., 511. 
35 [0039] Now, the convolution calculation, and its correspondence to Fourier transform performed in the filtering means 
12 is described, taking the case for extracting features of ridge line pitches as an example, 

[0040] When image data are expressed by g(x, y) (x = 0, 1, .... m - 1 and y= 0, 1 , n-1) and a filter is expressed 

by ft(x,y) (x= 0, 1 m h A and y=0,1 ,.2,..., n h - 1), the convolution z(x, y) of the image data sr(x,y) and a filter tyx.y) is 

calculated as following equation (1). 



40 



45 



50 



*(*, v) = £ Z 9(x- h y - j) Mi J) W 

[0041] This convolution calculation on a real plane is equivalent with a product calculation represented by following 
equation (2) of each frequency component on a Fourier transformed plane. 

F{g(x, y)*h(x, y)} = F{g(x, y))*F{h(x, y)} 

= G(u,v)*H(u,v) (2) 



55 where, A*B means convolution of A and B, F{A} means two-dimensional Fourier transform of A, and G(u, v) and H(u, 
v) mean the Fourier transforms of <7(x,y) and h(x,y) expressed on the Fourier transformed plane respectively. 
[0042] The two-dimensional Fourier transform is calculated as following equation (3). 
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10 



15 



25 



30 



35 



F{g(x, y)} - (?(«, «) = iEE *<*. V) e*P (-i2*(2f + ^)) (3) 

vmn x= q y=o \ ran/ 



where l/=0, 1 m- 1 and v = 0, 1, .... n-1. 

[0043] And by following equation (4), inverse Fourier transform, a Fourier transform is transferred on the real plane. 



g( Xi y) = F- l {G(u,v)} = -JL Z "f) G(u, v) exp (j2ir(£f + ^))(4) 

v mn «=o t?=o \ m n / 



[0044] Here, for shifting the direct current component G(0,0) to a center of the Fourier transformed plane, Fourier 
transform of (-1)* + >'<7(x f y) , hereafter called an optical Fourier transform, is often used for convenience of calculation. 
[0045] Returning to FIG. 1 , preparation of filters for extracting features of ridge line pitches in the filter file 15 will be 
20 described. 

[0046] For preparing a filter which filters out images of pitch frequency under and over fy, an ideal band-pass 
filter H(u, v) represented by following equation (5) and illustrated in FIG. 2 is designed on an optical Fourier transformed 
plane. 



_ a when R x < (ti* + v*)W < R 2 

H(u,v) = { (5) 

0 otherwise. 

[0047] Here, a is a constant for normalizing outputs of different filters, determined as follows, for example. 

a = V(F? 2 - fl?) (6) 



[0048] When the resolution of an image is 20 pixels per mm, the ridge line pitches are in a range from about 5 pixels 
to 20 pixels. This range from about 5 pixels to 20 pixels corresponds to a frequency range from about 100 cycles to 
25 cycles, difining 512 pixels as a cycle. 

[0049] This frequency range from 25 cycles to 1 00 cycles (a frequency band of 75 cycles) is divided into five ranges 
40 of frequency band of 1 5 cycles each, for example, for designing five filters corresponding to each thereof as shown by 
to h 5 in FIG. 3. 

[0050] Then, from these filters H,{u, v) (/= 1, 2, .... 5) designed on the optical Fourier transformed plane, inverse 
Fourier transforms are calculated for obtaining filters hfa, y) (/= 1 , 2, 5) to be applied for convolution calculation 
on the real plane. In this inverse Fourier transform, higher frequency components are ignored for simplifying the con- 

45 volution calculation. Thus, filters on the real plane are obtained approximating corresponding ideal band-pass filters 
to, to h 5 of FIG. 3. In the example, filter size of these transformed filters is 64 x 64 pixels. FIG. 4 shows shading images 
of these five filters defined by ft 1 to h s of FIG. 3. The data of these filters shown in FIG. 4 are stored in the filter file 15. 
[0051] In the filtering means 12, a convolution calculation expressed by the equation (1) is performed. When each 
of the five filters is denoted by h/x, y) (/ = 1, 2, 5), filtered data obtained through filtering by a filter /7/(x,y) are 

so expressed by Py(x f y) as follows, 

9r( x >y)= 9(x,y)*h f (x,y) (7) 

55 Five sets of filtered data p^x, y) (/= 1, 2, .... 5) are delivered to the image intensity calculation means 13. 

[0052] In the image intensity calculation means 13, image intensities P{i,m,n) of predetermined subregions C(m t n) 
of each set of the filtered data pXx.y) are calculated as following equation (8), 
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P(*\m,n)= E 9i(*,v) 2 (8) 

(x f y)6C(m l n) 

5 

[0053] Here, in the embodiment, each of the subregions C(m, n) is defined as a square region consisting of 32 x 32 
pixels as shown in FIG. 5, having its center on coordinates (m, n) arranged by every 16 pixels in x and y direction of 
the plane of 512 x 512 pixels, as represented by a bullet in FIG. 5, of each set of filtered data. 
[0054] The image intensity calculation means 13 delivers the calculated results P(i, m, n) to the feature extracting 
10 means 14. The feature extracting means 14 determine features of the predetermined subregions of the skin pattern 
image. When a filter h 1 (x.y), for example, represented by in FIG. 3 gives the maximum value of the image intensities 

P(i,m,n) (/= 1, 2 5) for a subregion C(m, n), ridge line pitch Qfeature t of the subregion C(m, n) is calculated as 

follows. 



(*i+* 2 )/2 
51? 

= ° = 1 5.7pixels (9) 
(25 + 40)/2 K v ' 

20 

[0055] Thus, features for every subregion are determined in the feature extracting means 1 4. 
[0056] Now, a second embodiment of the invention is described, wherein filtering of skin pattern images are per- 
formed by product calculation of each frequency component of the Fourier transforms. 

[0057] FIG. 6 is a block diagram of another filtering means 24 of the second embodiment corresponding to the filtering 
25 means 12 of FIG. 1 , comprising Fourier transformation means 21 for performing two-dimensional optical Fourier trans- 
form of image data read out from the image memory 1 1 , masking means 22 for performing product calculation between 
image data deliverd from the Fourier transformation means 21 and fitters prepared in a filter file 25, and inverse Fourier 
transformation means 23 for performing inverse optical Fourier transform of the masked data obtained from the masking 
means 22. 

30 [0058] The embodiment is described following the example applied in connection with the first embodiment for ex- 
tracting features of ridge line pitches of a skin pattern image. 

[0059] In the filterfile 25, the ideal band-pass filters H^u, v) (/= 1,2, ...5) designed on the optical Fourier transformed 
plane having data of 512 x 512 pixels are stored intact without transformation. The Fourier transformation means 21 
transforms the skin pattern image p(x,y) according to the equation (3) for obtaining an optical Fourier transform G(u, 
35 v) thereof. 

[0060] In the masking means 22, real and imaginary part of each frequency component of the Fourier transform G 
(u, v) are multiplied by those of the five band-pass filters Hi(u, v) one by one. 

[0061] Five sets of masked data G/u, v) = G(u, V) • H/u, v) are processed with inverse optical Fourier transform 
described in connection with the equation (4) in the inverse Fourier transformation means 23 for obtaining five sets of 
*o filtered data g&x, y) to be processed in the image intensity calculation means 13 of FIG. 1 as same as in the first 
embodiment. 

[0062] Filtering by convolution calculation as performed in the first embodiment needs a number of operations pro- 
portional to the filter size. So, the filters for convolution calculation are generally restrained in a size, 64 x 64 in the 
example described, accompanied with inevitable approximation. 

45 [0063] On the other hand, for filering by masking calculation as performed in the second embodiment, ideal band- 
pass filters can be desigend in a finite Fourier transformed plane with no approximation. Further more, for image data 
with its sides expressed by a power of 2, fast Fourier transform can be applied for high speed calculation. 
[0064] It is apparent that the optical Fourier transform can be replaced with the fast Fourier transform or an ordinary 
digital Fourier transform by preparing appropriate filters. Therefore, it can be said that there are more merits in the 

so second embodiment compared to the first embodiment when filter size is not sufficiently small. 

[0065] Now, a third embodiment will be described, wherein the feature extracting means 14 of FIG. 1 is replaced 
with another feature extracting means 30 illustrated in FIG. 7. 

[0066] In the third embodiment, there is defined an evaluation function, variable of which are feature parameters of 
each subregion, representing sum total of image intensities (multiplied by -1) of ail subregion and sum total of square 
55 sums of differences between values of the feature parameter of each subregion and that of neighbouring subregions. 
Each of the feature parameters which gives the minimum value of the evaluation function is determined as the feature 
of each subregion in the third embodiment. 

[0067] For the purpose, the feature extracting means 30 of FIG. 7 comprises minimizing means 31 and initial value 
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calculation means 32. 

[0068] From the image intensity calculation means 13 of FIG. 1 , the image intensities F\i, m, n) are delivered to the 
initial value calculation means 31 for determining initial values Sg(/n,n) of feature parameters, namely ridge line pitches 
that give a highest value of the image intensities fl(/,m,n) of each subregion C(m t n) as same as in the first embodiment. 
5 [0069] The map of the initial values s 0 (m,n) for all the subregions shows an initial two-dimensional feature pattern 
of the skin pattern image. For smoothing the initial two-dimensional feature pattern, the initial value So(m,n) of a con- 
cerning subregion C(m,n) is replaced with a most appropriate value s(m t n), among five values of feature parameters, 
which gives a minimum value of the evaluation function E(s{m,n)) represented by a following equation (10) in the 
minimizing means 31. 

10 

E(s(m y n)) = - £ P(s(m,n),m,n) 

(m,n) 

+ £ £ D(s(m,n) } s 0 {m\n')) (10) 
(m,n) (m' ( n') 

where, (m\n) means coordinates of neighbouring subregions of the concerning subregion C(m, nj,5x5*1 subregions 
20 shown in FIG. 8 in an example, and D(s(m, n), s^m', n)) is a function representing singularity of a value s(m, n), namely 
the current feature parameter, among its neighbourhood, expressed as follows in an example, a being a coefficient. 



D(s(m, n), s 0 (m', n)) = a{s(m, n) - s 0 (m', n*)} 2 (1 1) 

25 

[0070] Thus, the minimization of the evaluation function E(s(m t n)) is performed for every subregion C(m, n) replacing 
initial value s 0 (m, n) by an appropriate value s(m,n) one by one, and these processes are repeated several times 
determined for converging the minimization as shown by a flowchart illustrated in FIG. 9. 

[0071 ] In step 41 , initial value of the evaluation function Eis calculated according to the equation (1 0) with the initial 
30 values So(m,n) of feature parameters. Then, in step 42, value of the evaluation function Eis recalculated for a subregion 
C(m, n) with other values s M0 (m, n) one by one (checked in step 45) for replacing the initial value So(m, n) (in step 44) 
when the evaluation function Eis found reduced (in step 43). Steps 42 to 46 are performed for all subregions C{m, n) 
through step 46 and repeated several times through step 47. 

[0072] Thus a smoothed two-dimensional feature pattern giving the minimum value of the evaluation function E is 
35 obtained in the embodiment. 

[0073] Now referring to FIG. 1 0, a fourth embodiment providing a skin pattern image processor is described, wherein 
the extracting method of skin pattern features of the invention is applied. 

[0074] In the block diagram of FIG. 10, pixel value calculation means 75 are comprised in addition to the block 
diagram referred to in connection with the first, second or the third embodiment, for generating image data of a feature 
40 enhanced skin pattern. 

[0075] Taking the example of FIG. 5, each pixel belongs to four subregions since the subregions are defined over- 
lapped. Therefore, at most four values, are obtained for each pixel from filtered data sets corresponding to features 
assigned to the four subregions where the pixel belongs. Average of these at most four values for each pixel is output 
as each pixel value of the image data of the feature enhanced skin pattern image. 

45 [0076] In a fifth embodiment illustrated in FIG. 11 , skin pattern image processors of FIG. 10 are cascaded for en- 
hancing plural kinds of features of skin pattern image data. An example of the embodiment of FIG. 11 consists of two 
skin pattern image processors 80 and 81 , wherein skin pattern image data enhanced of ridge line pitch features by the 
skin pattern processor 80 are stored in a second image memory 806 to be processed with filters prepared for extracting 
ridge line direction features of the skin pattern image data. 

50 [0077] Therefore, a more smoothed skin pattern image can be obtained with features of two kinds enhanced in the 
embodiment. 

[0078] It goes without saying that feature data of each kind, to be used for classifying skin pattern images for example, 
can be obtained also in the embodiment from each feature extracting means 804 and 809 in addition to the feature 
enhanced skin pattern. 

55 [0079] The embodiments are described heretofore in connection with the example provided of five filters for extracting 
ridge line pitch features. However, any other appropriate filters can be used for extracting desired features of skin 
pattern images and processing them for enhancing the desired features thereof. 
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[0080] In the following paragraphs, other examples of filter sets prepared for masking Fourier transforms of skin 
pattern image data are described. 

[0081] FIG. 12 shows shading images of a filter set consisting of ten band-pass filters for extracting ridge line pitch 
features by masking optical Fourier transforms of skin pattern image data prepared in a first filter file of 811 of the 
5 embodiment of FIG. 11. Pixel values of each of the band-pass fitters of FIG. 12 are represented by following 
equation (12) symmetric with respect to its origin. 

" ^m)=r-(-^)//r-(-^)^ 12 ) 



15 where vector /-denotes coordinates of a pixel on an optical Fourier transformed plane of 512 x 512 pixels with origin 
at its center, r 0 and a r corresponding to a center frequency and a frequency width of pass-band respectively of each 
of the ten band-pass filters. 

[0082] In the examples of FIG. 12, r 0 are arranged by every 5.7 pixels from 28.4 to 79.6 pixels on condition o f is 1 2.0. 
[0083] By the way, shading images in FIG. 12 show center regions of the band-pass filters and the brighter dot 
20 represents the larger pixel value, darkest dot representing pixel value of 0. 

[0084] FIG. 13 shows shading images of another filter set consisting of 16 directional filters for extracting ridge line 
direction features by masking Fourier transforms of skin pattern image data prepared in a second filter file 81 2 of the 
example of FIG. 1 1 . Pixel values W e (i) of each of the directional filters of FIG. 1 3 are represented by following equation 
(13) line symmetric. 

25 

Wr) = exp(-^-«)') (13) 

30 

where, e e = (cos0; sinG) denotes unit vector of each direction 6 represented by each of the directional filters and a G , 
0.4 in the examples of FIG. 13, is a parameter representing phase angle width of the pass-band. 
[0085] As for the singulality function Dfsfm.n), So(m',rj)) in the evaluation function E of the equation (1 0), Euclidean 
35 distances expressed by following equation (14) are better to be applied here, for smoothing directional features by 
directional filtering, instead of the equation (11). 

D(s(m, n), s 0 (m', n)) = ot{(sin 29(r77, n) - sin 2Q(m\n f )f 

40 

+(cos 26(m, n) - cos 26(/n', n 1 )) } (14) 



where Q(m, n) is directional feature of a subregion C(m, n). 

[0086] Thus, by smoothing ridge direction features of a skin pattern image after smooting its ridge line pitch features, 
45 or vice versa, a more enhancing and reliable image processing can be provided in the present invention. 

[0087] Returning to FIG. 11 , two skin pattern image processors of FIG. 1 0 are cascaded in the fifth embodiment for 

obtaining smoothed skin pattern images with their features of two kinds enhanced or extracting two kinds of features. 

However, by preparing filters according to kinds of features to be extracted or enhanced, plural kinds of features can 

be extracted or enhanced of a skin pattern image with one skin pattern image processor. 
50 [0088] In a sixth embodiment having the same configuration with the fourth embodiment of FIG. 10, a filter set for 

extracting and enhancing plural kinds of features, ridge line pitch features and ridge line direction features for example, 

is prepared in the filter file 25. 

[0089] For preparing the filter set corresponding to the filter sets of FIG. 12 and FIG. 13, each filter is given from a 
product W Q ,(/) of a filter of FIG. 12 expressed by the equation (12) and a filter of FIG. 13 expressed by the equation 
55 (13) as follows. 
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W 9tr (r) - W 9 (r)W r (\r\) 



5 




10 [0090] If the ridge line pitch features are to be classified into ten values of r 0 and the ridge line direction feachures 
are to be classifyed into 16 values of 6 as same as FIGs. 12 and 13, a set of 160 filters is thus prepared. 
[0091] In the image intensity calculation means 13 of the sixth embodiment, image intensities P(6, r 0 , m, n) of the 
predetermined subregions C(m, n) are calculated for every 1 60 filtered data g[Q, r 0 , x, y) obtained from the skin pattern 
image <7(x,y) by filtering with the 1 60 filters W Q r (!) as follows in a same way as the equation (8). 

15 

P(0,r 09 m y n)= £ 9{0,r 0} x,y) 2 

20 

[0092] In the feature extracting means 30 of the sixth embodiment, vector values 

representing ridge line pitch featuers and ridge line direction features for the predetermined subregions C{m, n) are 
obtained from the image intensities P(e, r 0 , m, n) calculated by the image intencity calculation means 13 in a similar 
way as described in connection with FIGs. 7 to 9. 
30 [0093] In the initial value calculation means 32, initial values v 0 (/77,n) are calculated from parameters (G, r 0 ) of filters 
W G J^i) which give maximum intencities P(G, r 0 , m,ri) of filtered data for the predetermined subregions C(m, n). Starting 
from the initial values v 0 (m,n), a two-dimensional feature pattern of the vector values v(m t n) is determined in the min- 
imizing means 31 to minimize an evaluation function represented by a following equation (16). 

E(v(m,n)) m - £ P(6 t r 0 ,m t n) 

(m,n) 

+ £ £ D(v(m t n) 7 vo(rn' t n')) 

40 (nfjnf) 

where, 

45 D(v(m t n),v 0 (m , ,n')) = a{(sin 2G(m,n) - sin 2Q(m' t n')) 2 

+(cos 26(m, n) - cos 26(m', ri)) 2 } 
+p{r 0 (m, n)- r 0 (m', n*)) 2 

50 

[0094] Thus, two kinds of features of each subregion of a skin pattern image are extracted in the sixth embodiment. 
The smoothed skin pattern image is obtained by the pixel value calculation means 75 in the same way with the fourth 
embodiment from filtered data corresponding to the two kinds of features extracted for each subregion. 
[0095] And as beforehand described; still more kinds of features can be extracted and still more kinds of filtering or 
55 smoothing can be performed with still more kinds of filter sets in the scope of the invention. 
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Claims 

1. An apparatus for extracting features of predetermined subregions of a skin pattern image, comprising: 

5 an image memory (11) for storing image data of the skin pattern image; 

a filter file (15) for storing a plurality of filters consisting of two-dimensional data prepared according to the 
features to be extracted; 

filtering means (12) for outputting filtered data sets, each of said filtered data sets obtained from said image 
data by convolution calculation with each of said plurality of filters; 
10 image intensity calculation means (1 3) for calculating image intensity of each of the predetermined subregions 

of each of said filtered data sets; and 

feature extracting means (1 4) for determining the features of the predetermined subregions of the skin pattern 
image referring to said image intensity of each of the predetermined subregions of each of said filtered data 
sets, wherein 

15 said feature extracting means (1 4) determine the features of the predetermined subregions of the skin pattern 

image according to feature parameters minimizing an evaluation function, said evaluation function being a 
linear function of said image intensity of each of the subregions of one of said filtered data sets corresponding 
to a feature parameter assigined to said each of the subregions, and a value corresponding to singularity of 
said feature parameter assigned to said each of the predetermined subregions compared to the feature pa- 

20 rameters assigned to certain of the predetermined subregions neighbouring to said each of the predetermined 

subregions. 

2. An apparatus for extracting features of predetermined subregions of a skin partem image, comprising: 

25 an image memory (11 ) for storing image data of the skin pattern image; 

a filter file (15) for storing a plurality of filters consisting of two-dimensional data prepared according to the 
features to be extracted; 

filtering means (24) for outputting filtered data sets, each of said filtered data sets being an inverse Fourier 
transform of a masked data set obtained from a Fourier transform of said image data by product calculation 
30 of each frequency component of said Fourier transform and each of said plurality of filters; 

image intensity calculation means (1 3) for calculating image intensity of each of the predetermined subregions 
of each of said filtered data sets; and 

feature extracting means (1 4) for determining the features of the predetermined subregions of the skin pattern 
image referring to said image intensity of each of the predetermined subregions of each of said filtered data 

35 sets, wherein 

said feature extracting means (1 4) determine the features of the predetermined subregions of the skin pattern 
image according to feature parameters minimizing an evaluation function, said evaluation function being a 
linear function of said image intensity of each of the subregions of one of said filtered data sets corresponding 
to a feature parameter assigined to said each of the subregions, and a value corresponding to singularity of 

40 said feature parameter assigned to said each of the predetermined subregions compared to the feature pa- 

rameters assigned to certain of the predetermined subregions neighbouring to said each of the predetermined 
subregions. 

3. An apparatus according to claim 1 or 2, further comprising pixel value calculating means for calculating a value of 
45 each pixel of the filtered skin pattern image from values of pixels of those of said filtered data sets that correspond 

to each of the features extracted of the predetermined subregions where said each pixel belongs. 



PatentansprQche 

50 

1. Vorrichtung zum Extrahieren von Merkmalen vorgegebener Unterbereiche eines Hautmusterbildes, mit: 
einem Bildspeicher (11) zum Speichern von Bilddaten des Hautmusterbildes; 

einer Filterdatei (15) zum Speichern mehrerer Filter, die aus zweidimensionalen Daten bestehen, die in Uber- 
55 einstimmung mit den zu extrahierenden Merkmalen vorbereitet worden sind; 

Fitterungseinrichtungen (12) zum AusgebengefilterterDatenmengen, wobeijededergefilterten Datenmengen 

aus den Bilddaten durch eine Faltungsberechnung mit jedem dermehreren Filter erhalten wird; 

einer Bildintensitat-Berechnungseinrichtung (13) zum Berechnen der Bildintensitat jedes der vorgegebenen 
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Unterbereiche jeder der gefilterten Datenmengen; und 

einer Merkmalsextraktionseinrichtung (14) zum Bestimmen der Merkmale der vorgegebenen Unterbereiche 
des Hautmusterbildes, die auf die Bildintensitat jedes der vorgegebenen Unterbereiche jeder der gefitterten 
Datenmengen Bezug nimmt, wobei 

s die Merkmalsextraktionseinrichtung (14) die Merkmale der vorgegebenen Unterbereiche des Hautmusterbil- 

des in Ubereinstimmung mit Merkmalsparametern bestimmt, die eine Bewertungsfunktion, die eine lineare 
Funktion der Bildintensitat jedes der Unterbereiche einer der gefilterten Datenmengen, die einem jedem der 
Unterbereiche zugewiesenen Merkmalsparameter entspricht, sowie eines Wertes, der einer Singularity des 
jedem der vorgegebenen Unterbereiche zugewiesenen Merkmalsparameters entspricht, ist, im Vergleich zu 

10 den Merkmalsparametern, die bestimmten der vorgegebenen Unterbereiche zugewiesen sind, die dem vor- 

gegebenen Unterbereich benachbart sind, minimieren. 

2. Vorrichtung zum Extrahieren von Merkmalen vorgegebener Unterbereiche eines Hautmusterbildes, mit: 

15 einem Bildspeicher (11) zum Speichern von Bilddaten des Hautmusterbildes; 

einer Filterdatei (15) zum Speichern mehrerer Filter, die aus zweidimensionalen Daten bestehen, die in Uber- 
einstimmung mit den zu extrahierenden Merkmalen vorbereitet wurden; 

einer Filterungseinrichtung (24) zum Ausgeben gefilterter Datenmengen, wobei jede der gefilterten Daten- 
mengen eine inverse Fourier-Transformation einer maskierten Datenmenge ist, die aus einer Fourier-Trans- 
20 formation der Bilddaten durch Produktberechnung jeder Frequenzkomponente der Fourier-Transformation 

und jedes der mehreren Filter erhalten wird; 

einer Bildintensitat-Berechnungseinrichtung (13) zum Berechnen der Bildintensitat jedes der vorgegebenen 
Unterbereiche jeder der gefilterten Datenmengen; und 

einer Merkmalsextraktionseinrichtung (14) zum Bestimmen der Merkmale der vorgegebenen Unterbereiche 
25 des Hautmusterbildes, die auf die Bildintensitat jedes der vorgegebenen Unterbereiche jeder der gefilterten 

Datenmengen Bezug nimmt; wobei 

die Merkmalsextraktionseinrichtung (14) die Merkmale der vorgegebenen Unterbereiche des Hautmusterbil- 
des in Ubereinstimmung mit Merkmalsparametern bestimmt, die eine Bewertungsfunktion, die eine lineare 
Funktion der Bildintensitat jedes der Unterbereiche einer der gefilterten Datenmengen, die einem jedem der 
30 Unterbereiche zugewiesenen Merkmalsparameter entspricht, sowie eines Wertes, der einer Singularity des 

jedem der vorgegebenen Unterbereiche zugewiesenen Merkmalsparameters entspricht, ist, im Vergleich zu 
den Merkmalsparametern, die bestimmten der vorgegebenen Unterbereiche zugewiesen sind, die dem vor- 
gegebenen Unterbereich benachbart sind, minimieren. 

35 3. Vorrichtung nach Anspruch 1 oder 2, ferner mit einer Pixelwert-Berechnungseinrichtung zum Berechnen eines 
Wertes jedes Pixels des gefilterten Hautmusterbildes aus Werten von Pixeln jener Datenmenge der jedem der 
Merkmale entsprechenden gefilterten Datenmengen, die aus den vorgegebenen Unterbereichen, zu denen jedes 
Pixel gehort, extrahiert worden ist. 

40 

Revendications 

1. Dispositif pour extraire des caracteristiques de sous-regions predeterminees d'une image de motif de peau, 
comportant : 

45 

une memoire d'image (11) pour memoriser des donnees d'image de I'image de motif de peau, 
un fichier de filtres (15) pour memoriser une piuralite de filtres constitues de donnees bidimensionnelles pre- 
pares en fonction des caracteristiques k extraire, 

des moyens de filtrage (12) pour delivrer en sortie des ensembles de donnees filtrees, chacun desdits en- 
50 sembles de donnees filtrees etant obtenu k partir desdites donnees d'image par un calcul de convolution & 

Paide de chaque filtre de ladite piuralite de filtres, 

des moyens de calcul d'intensite d'image (13) pour calculer I'intensite d'image de chacune des sous-regions 
predeterminees de chacun desdits ensembles de donnees filtrees, 

des moyens d'extraction de caracteristiques (14) pour determiner les caracteristiques des sous-regions pre- 
55 determinees de I'image de motif de peau se rapportant k ladite intensite d'image de chacune des sous-regions 

predeterminees de chaque ensemble desdits ensembles de donnees filtrees, caracterise en ce que 
lesdits moyens d'extraction de caracteristiques (1 4) determinent les caracteristiques des sous-regions prede- 
terminees de I'image de motif de peau en fonction de parametres de caracteristiques minimisant une fonction 
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devaluation, ladite fonction devaluation etant une fonction lineaire de ladite intensite d'image de chacune des 
sous-regions d'un ensemble parmi lesdits ensembles de donnees filtrees correspondant a un parametre de 
caracteristique attribue a chacune desdites sous-regions, et d'une valeur correspondant a une singularity dudit 
parametre de caracteristique attribue a chacune desdites sous-regions predeterminees compare aux para- 
metres de caracteristiques attribu6s a certaines des sous-regions predetermines avoisinant chacune desdi- 
tes sous-regions predeterminees. 

Dispositif pour extraire des caracteristiques de sous-regions predeterminees d'une image de motif de peau, 
comportant : 

une m£moire d'image {11 ) pour m6moriser des donnees d'image de I'image de motif de peau, 
un fichier de filtres (15) pour memoriser une plurality de filtres constitues de donnees bidimensionnelles pre- 
pares en fonction des caracteristiques a extraire, 

des moyens de filtrage (24) pour deiivrer en sortie des ensembles de donnees filtrees, chacun desdits en- 
sembles de donnees filtrees etant une transform6e de Fourier inverse d'un ensemble de donn6es masquees 
obtenu a partir d'une transformee de Fourier desdites donnees d'image en calculant le produit de chaque 
composante frequentielle de ladite transformee de Fourier et de chaque filtre de ladite plurality de filtres, 
des moyens de calcul d'intensite d'image (13) pour calculer i'intensite d'image de chacune des sous-regions 
predeterminees de chaque ensemble desdits ensembles de donnees filtrees, et 

des moyens d'extraction de caracteristiques (14) pour determiner les caracteristiques des sous-regions pre- 
determinees de i'image de motif de peau se rapportant a ladite intensite d'image de chacune des sous-regions 
predeterminees de chaque ensemble desdits ensembles de donnees filtrees, caracterise en ce que 
lesdits moyens d'extraction de caracteristiques (14) determinent les caracteristiques des sous-regions prede- 
terminees de I'image de motif de peau en fonction de parametres de caracteristiques minimisant une fonction 
devaluation, ladite fonction devaluation etant une fonction lineaire de ladite intensite d'image de chacune des 
sous-regions d'un ensemble parmi lesdits ensembles de donnees filtrees correspondant a un parametre de 
caracteristique attribue a chacune desdites sous-regions, et d'une valeur correspondant a une singularity dudit 
parametre de caracteristique attribue a chacune desdites sous-regions predeterminees comparee aux para- 
metres de caracteristiques attribues a certaines sous-regions predeterminees avoisinant chacune desdites 
sous-regions predeterminees. 

Dispositif selon la revendication 1 ou 2, comportant en outre des moyens de calcul de valeur de pixel pour calculer 
une valeur de chaque pixel de I'image de motif de peau filtree a partir de valeurs de pixels de ces ensembles parmi 
lesdits ensembles de donnees filtrees qui correspondent a chacune des caracteristiques extraites des sous-r6gions 
predeterminees auxquelles appartient chaque pixel. 
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FIG. 9 
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FIG. 12 
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FIG. 13 
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FIG. 14 PRIOR ART 
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